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T— H It is believed that soon after the Planck time, Einstein's general relativity theory should 

^ be corrected to an effective quadratic theory. Numerical solutions for the anisotropic 

^0 generalization of the Friedmann "flat" model _E 3 for this effective gravity are given. 

QQ It must be emphasized that although numeric, these solutions are exact in the sense 

QQ that they depend only on the precision of the machine. The solutions are identified 

^O asymptotically in a certain sense. It is found solutions which asymptote de Sitter space, 

Minkowski space and a singularity. This work is a generalization for non diagonal spatial 

^ _ ' metrics of a previous result obtained by one of us and a collaborator for Bianchi / spaces. 

CN Keywords: effective gravity; homogenous solutions; isotropisation 

1. Introduction 



The semi-classical theory considers the back reaction of quantum fields in a classical 
geometric background. It began about forty years ago with de WittP, and since then, 
its consequences and applications are still under research, see for example^. 

Different from the usual Einstein-Hilbert action, the one loop effective gravita- 
tional action surmounts to quadratic theories in curvature, see for example Refs. [TJ 
[3 HI It is the gravitational version of the Heisenberg-Euler electromagnetism. As it 
is well known, vacuum polarization introduces non linear corrections into Maxwell 
electrodynamics, see Schwingei 5 , first obtained by Heisenberg and Eulei 6 . 

This quadratic gravity was studied in the pioneering works of GEO and I 9 ! 10 ! 11 ! 
Starobinsky idea was that the higher order terms describe a slowly decaying dark 
energy i.e. to produce slow-roll inflation. It is of interest for example in the context of 
the final stages of evaporation of black holes, inflationary theories^, in the approach 
to the singularity^ and also in a more theoretical context!^. 
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The effective gravity was apparently first investigated by Tomita ^ for gen- 
eral Bianchi / spaces. They found that the presence of anisotropy contributes to 
the formation of the singularity. In Ref. [15] it is shown that the equations for the 
anisotropic part of the Bianchi type I metric can be integrated for a general f(R) 
theory including the quadratic R + R 2 gravity. Ref . fT6l shows that a quadratic Weyl 
theory is less stable than a quadratic Riemann scalar R 2 . Homogenous solutions in 
the context of quadratic gravity was also addressed b}l 18 l 19 l 2l ^lHi SchmidtP^l does a 
very interesting review of higher order gravity theories in connection to cosmology. 

We have previously numerically investigated Bianchi / type solutions^ for the 
diagonal line element only. Soon after this result we also investigated the stability 
of the this particular Bianchi / case^, also for the diagonal line element. We must 
mention that the question of stability of the de Sitter solution in quadratic gravity 
was previously addressed by ^51. It turns out that for zero cosmological constant 
A = 0, Minkowski geometry is structurally stable in the sense that there is a basin 
of attraction to Minkowski solution. And for A > de Sitter type geometry is 
structurally stable, also in the same sense. Thus, soon after the Planck era, the 
effective one loop quadratic gravity predicts for the very particular Bianchi I models, 
that there is a basin of attraction to Minkowski space for A = and a basin of 
attraction to de Sitter space for A > 0, see^. 

The purpose of this work is to consider more general geometries than the ones 
considered irPSancpS for instance with a non diagonal line element. Bianchi / is the 
anisotropic generalization of the Euclidean E 3 Friedmann "flat" space. The time 
like vector is orthogonal to the group orbit, geodesic and with zero vorticity. For an 
excelent and recent review in Bianchi cosmologies, see Ref. l2lfl 

For any spatially homogenous space-time, the dynamical equations of motion, of 
any metric theory, result in a non linear system of ordinary differential equations. 
For the particular quadratic gravity, the ordinary differential equation is of degree 
4. 

The initial conditions chosen are very near to exact solutions of the quadratic 
theory, and also only classical vacuum source is considered, which we believe is 
the best description soon after the Planck era. The solutions are obtained numer- 
ically and they are understood asymptotically in a certain sense described in the 
text. It is shown that depending on the initial conditions and parameters, the de 
Sitter solution, or Minkowski solution, or a singularity type solution is obtained 
asymptotically. 

The paper is organized as follows. Section 2 presents the field equations for the 
Bianchi / case, and also the Newman-Penrose coefficients. In Sec 3, three different 
initial conditions are analyzed: one that converges to de Sitter space, one that 
converges to Minkowski space, and one that converges to a singularity. We believe 
the index notation should be clear from the text. 
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2. Theory and Development 



The field equations for the semiclassical theory, are obtained performing metric 
variations in the gravitational Lagrangian 



C 



-A- /?-!- .. j R ab R ab - *i? 2 ] >/?- 



+ £ q , (i) 



where C q is the quantum part of Lagrangian and a and /3 are constants. For the 
spatially homogenous space they are described by the tensor E — E ab u a <E> u> b , 
LO A = dt, 

E ab = G ab + l -g ab K -U- l -c^j H { 2 - aH% = 0, (2) 

where 

Gab = Rab — ~^9abR, 

H a b = 7j9abR~ — 2RR a b — 2.g Q bV c V c i? + 2R. ab , 

H a b = y9abR C Red — V c V C i? a f, — -5afcV c V C i? 
+R;ab — 2i?° Rcbda- 

The counterterms R 2 , R a bR ab , A and R in M are precisely the ones necessary in 
order to obtain a finite vacuum expectation value of the energy momentum tensor, 
see for example Ref.1271 A theory without these counterterms is inconsistent from the 
point of view of the renormalization of the quantum field in C q . The renormalized 
vacuum expectation value of the energy momentum tensor is set to zero, which 
emphasizes the effects of a theory that should have been considered from the start. 
We are disregarding any classical contribution in this work. 

The finite contributions to the vacuum expectation value of the energy momen- 
tum tensor are known only for very particular situations, for example, split rank 
spaces, for which the heat kernel can be obtained exactly ^, then the point split- 
ting method gives (0|T o f,|0) exactly. These solutions are known as self consistent, 
Gab = (0|T a f,|0) and there are very few cases, for example Ref. [2U1 For another very 
interesting example of a self consistent anisotropic solution supported by vacuum 
polarization see Ref. |3D1 

As in any other metric theory, the covariant divergence of E ab must be zero 

V a E ab = 0. 

From V a E ab — it follows that if E i4: = and E ia = initially, then they 
will remain zero for any time, and act as constraints on the initial conditions. 
Consequently these constraints are checked to test the accuracy of the numerical 
results, while E a p = represent the real dynamical equations of the problem (see 
for example Ref.lBTl p. 165). 
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The homogeneous and anisotropic metric is supposed in the following way, 



/«i 2 W 



o 







9ab(t) = 



a 2 (t) a±{t) 



\ 



a 4 {t) a£(t) 
V -1/ 



(3) 



ds 2 = -dt 2 +g a p(t)dx a dxl 3 , 



(1) 



where g a p is the spatial part. The replacement of the above line element in (pi) results 
in a non-linear fourth-order ordinary differential equation system in the functions 
cn{t), z = 1,2,3,4: 



d 3 



(5) 



First let us emphasize that every Einstein space, satisfying R a b = Ag a b/2, is an 
exact solution of this effective theory given in |2]) . Note the particular case when the 
constant A = 0: Vacuum solutions of Einstein's equations are also exact solutions 
of pj . So there is the following exact solution of Q for Bianchi / with the diagonal 
metric g a p(t) = diag(a 1 2 (i),a 2 2 (t),a 3 2 (t)) 



oi(*) = CeVA/« 

a 2 (t) = Ce 1 ^ 1 
o 3 (t) = C e f ^^, 



(6) 



where C is an integration constant. The particular case when A — > 0, is of course 
Minkowski space. 

The intention is to characterize the obtained numerical solutions in some way. 
The Weyl tensor follows from the Riemann tensor as 



Cabcd = Rabcd ~ ^{9caRbd + 9dbRca ~ 9cbRda ~ 9daRcb) 



R{g C a9db — 9cb9da)- 



A complex null basis can be defined 

k k a = K t a = k t a = l l a = I t a 

l L a K l a 1, 

with the corresponding null metric 



t L a b b a i Z a U, 



a r>b 



9AB = 9abA a B 



/ -1 00\ 

-10 00 

1 

V o o 10/ 



(7) 



(8) 



(9) 
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where A a and B b are the null vectors in (J8J). The Newman-Penrose complex coeffi- 
cients are in fact the tetrad components of the Weyl tensor 

ipo = C abcd k a t k c t , 

4>i = c abcd k a i b k c t d , 
i>2 = c abcd k a t b i c i d , 

^3 = C abcd k a l b t c l d , 
lp4 = C abcd t a l t°l . 

Also the Ricci tensor can be projected to the null tetrad, resulting into the inde- 
pendent components, 

Rkk = Rabk a k b 
Rkl = Rabk a l 
Rll = Rab^l 

Rkt = R ab k a t 

Rlt = Rabl a t 
Rtt = Rabt a t 

Ru = R ab t a ?, (10) 

the first 3 and the last one are real, and the remaining ones are arbitrary complex 
numbers comprising of course 10 independent components. 

The Newman-Penrose coefficients are related to the Petrov classification as 
shown for example in Refs. EH and [SB When all the i/j's are zero it's a Petrov 
type O, which characterizes a conformal Minkowski space and the Weyl tensor van- 
ishes. An exactly conformally flat solution is either a generalized Friedmann solution 
or an interior Schwarzschild solution as it can be seen for example in Ref. [32] pg. 
413. Minkowski and de Sitter geometries, are particular cases of conformally flat 
solutions. 

Before getting into the numeric result, we call the attention to the boolP^I, 
pp. 62- 64. As it is well written there, there is not a complete statement as to 
what constitutes a minimal set for ensuring that a cosmological model is close to 
Friedmann-Lemaitre model. This book is based on an article by Stoeger et al. ^ 
for which some assumptions are made: i) Einstein's equations are satisfied, ii) the 
source is a mixture of radiation and dust. 

We are not concerned if the two above conditions are satisfied in this present 
work. So we do not expect that their result, albeit being very interesting, should be 
verified in the particular context being discussed here. 

In the following we will present a particular example. Consider the Bianchi / 
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metric written in the appropriate coordinate base dx, dy, dz, dt 



9ab 



fa(t) 2 \ 
b(t) 2 
b(t) 2 

V -1/ 



(11) 



The time like vector u a = (0, 0, 0, 1) is geodesic and orthogonal to the group orbit. 
The magnetic part E ab = 0, and the electric E ab part of the Weyl tensor is 



e: 



/2iJj 2 0\ 

o -v> 2 o 

o o -v>2 o 

V o o oo/ 



where ip2 is the above defined Newman-Pcnrosc coefficient 

1 f- (b(t)) 2 +b(t)b(t) \ 
b(t) 2 J ' 



1p2 = 



() 



Also regarding (11), the expansion H = 1/36, = V c u c is 

H - 
Now consider 

a(t) 

b(t)-- 



1 a(t) 2 b(t) 
3aW + 3b(tj' 

A 
= l-t-i/3 s i n (t) 

A 

1 - £-V4 COs(i) ' 



(12) 



(13) 



(14) 



From ( 14 ) it is easily seen that ( 11 ) as t — > oo is 



lim g ab = 

t— ^oo 



/A 2 \ 

A 2 

A 2 

\ -l) 



(15) 



which is Minkowski geometry. 



Considering (|14|) and (|ll|), now we can explicitly obtain the asymptotic t — > oo 

2j) a 

i- 1/4 sin(t) 



limits of the Newman-Penrose coefficient ( 12 ) and the expansion ( 13 1 to be 

1 
C 



lim ip2 

t— >-oo 



lim E ab E ab 



h-^smtt) 2 
6 v ; 



lim H = -r 1/A cos{t) 



(16) 
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This example shows that Definition 2.1 of^ 



VH ab H" b /H 2 « e 
^E ab E ab /H 2 « e 



(17) 



< e << 1 is not necessary for the asymptotic approach to a Friedmann-Lemaitre 
model. Of course the authors of Ref. [33] were aware of this fact as it can be read on 
pg. 64 where they carefully write that: if Definition 2.1 is satisfied, then the metric 
can be locally written in an almost-RW form. 



As a final word we will use a different criteria, namely, inspired in (16), the 
asymptotic tetrad components of the Weyl and Ricci tensor 



lim C abcd « e 

t— »-oo 

Km R a b — > const. 

t— >oo 



(18) 



where lei << 1, and the const, values do not have to vanish, which is less restrictive 



than (17 1. We emphasize that the asymptotic behavior of the tetrad components 



(18) is a weaker criteria than the one given in Definition 2.1 of Ref. 33 

For instance, Kasner's^ exact solution and the exact Bianchi VII a gravity wave 



of Ref. 25] obey (18) with R ab = 0: on the other hand these two exact solutions do 



not satisfy (17). So in principle our criteria would not distinguish between Kasner 



exact anisotropic solution, and for example, the limit as t — > oo given in (15), with 



( 14 ) and ( 11 ), which is clearly a space for which isotropisation occurs in the strong 



sense. 



3. Numerical Solutions 

3.1. Asymptotically de Sitter Solution 

We choose the an initial condition near the exact de Sitter (pi) solution with t = 
and C = I, a — 1, /3 — —5.0 and A = 0.06. The only non null coefficients consistent 
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with the -E44 = constraint are 

01 (t) =(7 + 0.1 
oa(t) = C 

03 (t) = <7 

04 (t) = 0.2 
a"i(t) = CV^V6 

02 (t) = C V / A76 
a 3 (t) = Cv/A/S 
di(<) = CA/6 
d 2 (t) = CA/6 

03 (t) = CA/6 
a 1 (i) = C(A/6) 3 / 2 
a 2 (t) = C(A/6) 3 / 2 

a 3 (i) = 8.31792283 x 10~ 4 (19) 

The time evolution of the non zero Newman-Penrose coefficients, V'o-V'4 are 
shown in Fig. [I] For the line element chosen ipi = ^3 = 0. It can be seen that 
asymptotically the Newman-Penrose coefficients all vanish so the Weyl tensor is 
asymptotically zero. 

The time evolution of the non zero components of the Ricci tensor according 
to the null base Qare shown in Fig. [2j For the particular metric chosen, p| -Q, 
Rkt — Rit = 0. It is also shown in Fig. [2] that asymptotically, the Ricci tensor is 
a constant proportional to the null metric Q. The only consistent result with the 
field equation p| is 
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Fig. 1. Numerical evolution of the real and imaginary parts of the Newman-Pcnrose coefficients 
•ipO - ip4- For the initial condition given in the text. The numerical integration was done until the 
proper time t = 300. 
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Fig. 2. Numerical evolution of the non zero real and imaginary parts of the null tetrad J81 
components of the Ricci tensor \W\ . The numerical integration was done until the proper time 
t = 300. It can be seen that the Ricci tensor is asymptotically proportional to the null metric d9l. 



Wc have numerically checked up to t = 500 in proper time that R a t, = A/2g a i, 
with one part in 10 12 . Also the constraint -E44 = 0, is numerically verified in Fig. [3] 
which is a strong indication that the numerical result should be trusted. 



March 1, 2013 23:25 WSPC/INSTRUCTION FILE article 



Bianchi I solutions of effective gravity 11 



-I 




Fig. 3. Numerical evolution of the constraint -E44, indicating that the numerical solution is accu- 
rate. The integration was done for the initial condition and parameters specified on the text and 
until the proper time 300. 

3.2. Asymptotically Minkowski Solution 

In the numerical solutions the following values for the parameters C — 1.0, a = 2.0, 
/3 = —5.0, A = 0, where chosen, and initial conditions (t = 0) near Minkowski 
solution 

01 (t) = C + 2 
a 2 (t) = C 

03 (t) = C 

04 (t) = 2.0 (20) 

The only non null initial conditions for the derivatives are given in Table [l] and 
as before, the numerical value of 0*3 is consistent with the constraint £^4 = 0. 

Table 1. The only non null initial con- 
ditions together with d20t , consistent 
with the E44 constraint 



ai(0) 03(0) 



a 3 



1.0 



0.2 



-2.45214348 x 10" 



First of all, linearization of (12| on a slowly varying background geometry, S a f, << 
uj gives the following eigenvalues 



9ab = Vab ~ 

uj = 0, u> 



Sane 1 " 1 
1 



(21) 



in Riemann local coordinates where the components \S a b\ << 1. 

The time evolution of the non zero Newman-Penrosc coefficients, 1JJ0-1P4 are 
shown in Fig. [4] For the line element chosen ipi = ^3 = 0. The time evolution of the 
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Fig. 4. Numerical evolution of the real and imaginary parts of the Newman-Penrose coefficients 
ipo - Tp4- For the initial condition given in the text. The numerical integration was done until the 
proper time t = 200 . It can be seen directly from the graphs that u) = 2tt/T as 0.706, in good 
agreement with E^ w = l/y/a = 1/V2 w 0.707 



non zero components of the Ricci tensor according to the null base (fSj) are shown 
in Fig. [5} For the particular metric chosen, ^ -Q, R kt — R it = 0. 
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Fig. 5. Numerical evolution of the non zero real and imaginary parts of the null tetrad qSn 
components of the Ricci tensor JTM . The numerical integration was done until the proper time t = 
200. Again, form the graphs the higher frequency u) = 0.707. We have spe cifically checked that the 
slower frequency u) = 2tt/T pa 0.1820 is also in good agreement with J2l|, u> = 1/ \f— 6/3 ~ 0.1825. 



The time evolution of the constrain S44 is shown in Fig. [6] 
Since all the i\i — > 0, Fig.Hlthe Weyl tensor C% cd — > 0, and since the Ricci tensor 
R a b — > Fig. [5j then the Riemann tensor R a bcd — > 0. 

The numerical integrations shown in Fig. [4] and Fig. [5] were carried up to times 
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Fig. 6. Numerical evolution of the constraint E44, indicating that the numerical solution is accu- 
rate. The integration was done for the initial condition and parameters specified on the text and 
until the proper time t = 200. 



t = 1 X 10 6 showing that null tetrad components of the Riemann tensor R a bcd —> 
in one part in 1 x 10 6 and smaller values. It is in this sense that the solution is 
understood to asymptotically approach Minkowski space. 



3.3. Singularity 

Using exactly the same parameters of the preceding section C — 1.0, a 



2.0, 



P = —5.0, A = 0, a slightly different initial condition is chosen 

01 (t) = C + 2.0 

02 (t) = C 

03 (t) = C 

a 4 (i) = 2.0, (22) 

and the only non null initial conditions for the derivatives are given in Table [21 

Table 2. The only 
non null initial con- 
ditions together with 
1 ] '2 2 [ ) , consistent with 
the £44 constraint 



a 3 (0) 



a 3 



2.0 



-7.111111, 



This initial condition evolves very fast to a singularity characterized by the 
increase of the curvature scalars R a bR ab and R a bcdR abcd shown in Fig. [7] The con- 
straint which should be zero -E44 shown in Fig. [51 prove that the numerical result is 
accurate. Note that as the singularity is approached, the numerical errors increase, 
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Fig. 7. The increase of the scalar curvature invariants R a i,R ab , R a bcdR abcd characterizing a 
singularity. 

which is expectable. Also, the numerical integration can be carried further and fur- 
ther, and the constraints are not really satisfied at all, and is not shown since we 
believe the result should not be trusted. 
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Fig. 8. The constraint -E44. Note that as the singularity is approached, the numerical errors 
increase, which is expectable. 
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4. Conclusions 



In this present work, the numerical solutions of Bianchi / type are analyzed in the 
effective gravity context. This is the anisotropic generalization of the Euclidean E 3 
Fricdmann "flat" space. The quadratic gravity is the result of vacuum polarization 
counter terms which must be introduced into Einstein's theory of gravitation. It 
should be the most natural theory just after the Planck era. It must be emphasized 
that although numeric the solutions are exact in the sense that they depend only 
on the precision of the machine. 

We have previously, numerically investigated Bianchi / type solutions ^3. We 
found that there is a basin of attraction to Minkowski space for A = and a basin 
of attraction to de Sitter space for A > 0, see ^2. In this sense, Minkowski and 
de Sitter solutions are structurally stable according to the effective gravity, for the 
particular Bianchi I models we analyzed. 

The numerical solutions are given in section 3. In section 3.1 the initial condition 
is chosen near the de Sitter exact solution. In section 3.2. and 3.3 the the solution 
is chosen near Minkowski exact solution. 

The solutions are understood asymptotically in the following sense. The Rie- 
mann tensor depends on the Weyl tensor, Ricci tensor and Riemann scalar. The 
numerical solutions are characterized asymptotically according to the null tetrad 
components of the above tensors. 

In this sense, for the first solution under consideration in section 3.1, the Weyl 
tensor is zero and the Ricci tensor is proportional to the metric R a b = A/2g a i, 
asymptotically and we identify this solution with de Sitter space. In section 3.2 the 
solution asymptotes Minkowski space in sense that the Weyl tenosr and Ricci tensor 
vanish asymptotically. The convergence to Minkowski space is not so fast and we 
shall not address this question in this present work. In section 3.3 a solution that 
converges to a singularity is presented. 

In all the cases we considered, the numerical behavior described in the text was 
checked for much larger times that the ones plotted, and we believe that asymptot- 
ical interpretation is the correct one. 

So this work is in accordance with previous results obtained by on us for the more 
particular case of the diagonal Bianchi I caseEES. As in the diagonal metric case, 
note the presence of the singularity, depending on the initial condition: this theory 
certainly can not be understood as a complete one because an initially reasonable 
condition evolves to a singular universe. On the other hand, depending on the initial 
condition isotropisation occurs in a weak sense, see discussion on pg. 3 of this work. 

Since the initial conditions are near de Sitter and Minkowski solutions, we can 
speculate that also for Bianchi I solutions, de Sitter and Minkowski space should 
be structurally stable according to the effective gravity for the non diagonal metric 
also, in the sense that there should be basins of attraction to these solutions. We 
intend to address this stability issue in a future work. 
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